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We carried out a Pade approximant analysis on a compact factor of the T- matrix for NN scattering 
to explore the nonperturbative renormalization prescription in a universal manner. The utilities and 
virtues for this Pade analysis were discussed. 
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O . INTRODUCTION 

o ■ 

Since Weinberg's seminal work[l|, the effective field theory (EFT) approach to the nuclear forces has been exten- 
sively investigated [2|. However, such applications are plagued by severe nonperturbative UV divergence. For the 
EFT approach to be useful, the regularization and subtraction scheme must be carefully worked out together with 
a consistent set of power counting rules. In other words, appropriate renormalization pres crip tion is needed in non- 
perturbative regimes. There have been many contributions to this issued, EL S S H E EL 3 13, 14, 15, 16| 
(More could be found in Ref.Q), among which there are some controversies and debates. At some points, different 
Tj- ■ approaches could lead to rather disparate predictions [17]. 

Recently, a compact parametrization of the T-matrix is proposed in Ref.[18j|. with which the obstacle for renor- 
malization being identified as the compact form of the T-matrix. In a concrete example [l9j]. it was shown that the 
T-matrix could only be renormalized through 'endogenous' counter terms, which result in nontrivial prescription 
dependence. Such prescription dependence must be removed or fixed by imposing appropriate physical boundary con- 
t-H ' ditions, for instance, through certain procedure of data fitting, as is frequently done in literature. The conventional 
"^t" 1 power counting could be preserved within such procedures. 

In this short report, we sketch a Pade approximant analysis basing on the aforementioned par ametrization for the 
. T-matrix. The paper is organized as follows: In Sec. II, the parametrization proposed in Re.f[18l| is briefly described 
"^T* \ with some related remarks. In Sec. Ill, the Pade approximant of a factor in the compact parametrization of the 
T-matrix is employed to parametrize the nonperturbative prescription dependence. Then predictions for phase shifts 
are made at various chiral orders and Pade approximants. The instrumental utilities and other aspects of this analysis 
d ' will be discussed in Sec. IV. The report is summarized in Sec. IV. 
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THE COMPACT PARAMETRIZATION 

To describe low energy NN scattering, one first constructs the potential V from xPT[H up to certain chiral order, 
then computes the T-matrix through Lippmann-Schwinger equation (LSE), which is plagued with severe nonpertur- 
bative UV diver gen ces. To appreciate the crucial aspects, a compact form of T-matrix in diagonal channelsjlO] is 
proposed in Ref.[l8( as follows basing on LSE, 

7f J-, = 77 1 1 ; \ G2s+i Lj {p',p;E), (1) 

T2s+i Lj {p',p;E) V2s+i Lj [p',p) 

/kdk 2 
——V2s +lLj (p', k)G (k; E)T2 S+1Lj {k,p; E) 

V2s+i Lj {p',p)T2s+i Lj {p',p;E) 

Here G (k;E) = 1/(E - k 2 /M + ie), with M bein g nu cleon mass, p' and p being the off-shell external momenta. 
Using the on-shell relation between K- and T-matrix [2 lj: 1/T(p) — l/K(p) +ijr-p, we arrive at the following on-shell 
relations (from now on, we omit the subscript ' 2S+1 T,/') 

G( P ) =V- 1 (p)-K- 1 (p)-i^p; (3) 

T-\p) =V-\p)-G(p). (4) 

Obviously, Q assumes all the nonperturbative divergences in a compact form. Any approximation to the quantity Q 
leads to a nonperturbative scheme for T. Here we should remind that the power counting is applied in the construction 
of the potential V(p,p'). 
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In perturbation theory, UV divergences are removed order by order before the amplitudes are summed up. While 
for Eq.{T]) or (jT0|) in nonperturbative regime, infinitely many UV divergent amplitudes like VGnVGnVGn ■ ■ ■ V must 
be lump summed into a compact form. Then one must specify the order for implementing the following two in- 
commutable procedures: subtraction versus nonperturbative summation. So a natural discrimination arises between 
'endogenous' and 'exogenous' counter terms (or equivalent operations) that are introduced before and after the sum- 
mation respectively[18|, The compact form of T-matrix fails the 'exogenous' counter terms [l8L fl9| ] . In other 
words, the renormalization through 'endogenous' counter terms is the only sensible procedure in nonperturbative 
regime. Since the Schrodinger equation approach[llj is intrinsically nonperturbative, any successful subtraction in 
the Schrodinger equation approach serves as a concrete instance for 'endogenous' renormalization in action [11113]. In 
practice, 'endogenous' subtraction is a formidable task in the T-matrix formalism: To complete the subtraction AND 
summation to ALL orders! That is why various forms of finite cutoff prevail in literature. In whatever approaches, 
the compact form of T-matrix persists and makes nonperturbative prescription dependence strikingly different from 
the perturbative cases 1|| 22|: Physical boundary conditions must be imposed as a nontrivial procedure. 

In what follows, all the EFT couplings are collectively denoted by [C..J and pion mass by m n . In general, any 
prescription could be parametrized by a set of dimensionless constants [q...] and a dimensional scale including 
various finite cutoff schemes. 



A PADE- AIDED ANALYSIS OF THE EFT FOR NN SCATTERING 



Motivation 



From the above discussions, it is clear that the nonperturbative renormalization prescription is solely assumed in 
the factor Q. Obviously, the Q factor could not be perturbative in terms of the EFT couplings, and its nontrivial 
nonperturbative prescription de pen dence is parametrized by the parameters to be physically fixed. These points 
have been demonstrated in Ref.jl9(. From the point of view of Eq.((T]), the main issue in the EFT for low energy 
NN scattering is to work out an appropriate procedure or prescription of renormalization in the nonperturbative 
regime so that the EFT power counting schemes (encoded in V) remain intact. This requires the full solution of the 
nonperturbative factor Q. 

Now we need a general formalism to describe the nontrivial features of the renormalization in the nonperturbative 
regimes. Due to the difficulty in obtaining the full analytical nonperturbative solutions, certain approximation must 
be employed in practice, including various numerical approaches^. Here, we wish to propose an approximation 
approach that is, we feel, more analytical in the nonperturbative regime. The starting point is just the compact 
parametrization defined in Eq.(fT]). We should note that in this parametrization the direct EFT component is the 
potential, which is constructed using EFT power counting. Thus the EFT is naturally incorporated in the following 
analysis through the potential. We will return to this issue in the third subsection. 



Formulation 



The idea is very simple, we parametrize the factor Q in terms of Pade approximant. This is reasonable as Q is 
nonperturbative in terms of the EFT couplings and prescription parameters. In formulae, we employ the following 
parametrization of the Q factor (p = V ME) : 

nl v., f vo + v lP 2 + --- M . \ 2 M . 

S(p)||pade => \ s + s ip z +77. ~ J^ ip j " Taylor ^ - 9 (°) + 9 W p + J^ tp ' ( ' 

Here, the Taylor series is also listed as an expansion in much lower energy regions. Note the significant distinction 
between the Pade analysis of the Q factor here and that of the whole T-matrix: In the former case, EFT is indispensable 
in the systematic construction of the kernel (potential), while in the latter case, EFT plays no role at all. 

Obviously Q assumes all the prescription dependence through vi , Sj or that are nontrivial functions of the EFT 
couplings [C...,M, m w ], and the prescription parameters, [q...,fi\. Then instead of [(?..., /i], we could use i>i,6j or g< n -\ 
to parametrize the renormalization prescription of the T-matrix within the EFT approach, 

T-\p; [C..,M,m w ]; [i/.. .,*..]) =► V^ 1 (p; [C..,M,m„]) - ^ T 7^2 ^ ' " ' + 7~*P' ( 6 ) 

On + Olir + • ■ ■ 471" 
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T \ P ; [C...,M, m*]-, => V~ l {p\ [C...,M, mj) - (g (0) + g {1)P 2 + •••) + (7) 

Intuitively, nonperturbative prescription dependence in , <5j or <?(„) could be understood from the rigorous non- 
perturbative solution of on-shell T-matrix for So channel scattering with contact potential at next-to-leading order 
(Nlo) and next-to-next-to-leading order (Nnlo) [19(] : 

Nlo:V (2) = Co + C 2 (p 2 + p' 2 ); 

Nnlo : V (4) = C + C 2 {p 2 +p' 2 ) + C 4 p 2 p' 2 + C 4 (p 4 + p' 4 ); 

rp-i N ° + N lP 2 + N 2 p* M , 

= D + D lP i + D 2 p± + D^ + - h + -^ (9) 

with [Ni, Dj], which correspond to [z/j, Sj], being polynomials in terms of coupling [C...] and [</„]. Here the constants 
[J n ] come from divergent loop integrals and hence parametrize the renormalization prescription, like the set [q..., /i] 19]. 
Then Ni, Dj effectively parametrize a prescription due to their nontrivial dependence on [</„]. In Pade approximant, 
i>i, 5j or g(n) take over the role of N iy Dj, and the physical boundary conditions are to be imposed on z^, Sj or <?(„). 



Power counting and nonperturbative prescription 

As is stressed in the previous section, EFT and its power counting rules enter through the potential. Through 
Eq.([T]) the EFT elements and their power counting rules are carried over to the T-matrix in a nonperturbative manner. 
Since the kernel of the LSE or the potential is perturbative, its renormalization is still perturbatively implcmentable 
within the EFT power counting rules. Therefore, for the T-matrix, Eq.{T]) just provides a concise separation of 
the nonperturbative renormalization information from other things. Both Q and V are in principle EFT objects, 
the sole and important distinction is that the former exclusively carries the information about the nonperturbative 
renormalization prescription. Thus, we could append some subscripts to Eq.([T]) as follows, 

r EFT,n.p.t» =^EFT,p.t.b)-^EFT,n.p.t.(p)- (10) 

In general, within a natural EFT, the nonperturbative objects (e.g. Q) should also exhibit certain degree of 
naturalness in the sense that, the scales involved in such objects should not deviate very much from the natural sizes. 
But the renormalization in nonperturbative regimes does allow for other unconventional scenarios, without violating 
natural power counting rules [19(. As a matter of fact, the Pade analysis of Q alone according to Eq.{T]) does not 
affect anything of the EFT power counting rules, i.e., the implementation of the EFT power counting rules and the 
nonperturbative renormalization procedures are disentangled. The possible subtleness in the nonperturbative factor 
Q should not be misunderstood as the inconsistency or the unnaturalness of the EFT power counting rules or even as 
the inapplicability of EFT method at all. That is, the unnaturalness in i/j , Sj or could have nothing to do with 
the inconsistency of EFT power counting. 



Fitting and predictions: So channel 

With the preceding preparations, we can demonstrate the predictions of phase shifts for the 1 S'o channel NN 
scattering at different orders of potential using different Pade approximants. For each case, the prescription param- 
eters, i.e., the Pade parameters are fixed through fitting the phase shifts in the low energy ends. The laborious loop 
integrations and 'endogenous' subtractions are naturally avoided. 

We will employ the potentials that are worked out in Ref.Q (denoted as EGM from now on), which contains no 
energy dependence, and less contact couplings-a favorable aspect for fitting the 'physical' values for Vi,5j or <?(„). 
One could well employ other construction schemes for potentials. In fact, one could compare any pair of potential 
schemes only through fitting the parameters z/$ , Sj or . It is obvious that at any chiral order with any Pade 
approximant, different Pade parameters would yield rather different phase shifts curves. We will not show the figures 
for demonstrating such nontrivial prescription dependence due to space limitation. Let us focus on more interesting 
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figures with the Pade parameters determined through boundary conditions: fitting in the low energy regions: (1) At 
Lo, Ti ab G (0,3) in units of MeV;(2) At Nlo , while T lab G (0.2,13); (3) At Nnlo, T lab G (3,23). The phase shift is 
obtained from the following formula, 

5(p)= a rctan{-^f-i- T -Re(0(p))') 1. (11) 



4vr \V(p) 

The phase shifts predicted at Lo, Nlo and Nnlo are depicted in FigJT](a), (b), (c) respectively. At each order, three 
Pade approximants are shown respectively: (1) Re(0) w g/o) (dotted lines); (2) Ke(Q) w gt \ + g(i)P 2 (dashed lines); 
(3) Re((?) w (l + ^P 2 )/(f^ + ^P 2 ) (solid lines). From these diagrams, one could find that, at each order, the prediction 
improves as more Pade parameters are present, which is a natural tendency. One could also anticipate that with each 
Pade approximant, the predictions should also improve as higher order terms are present in the potential, which are 
responsible for the interactions at higher energy. The results are shown in FigO In each figure, the predictions are 
compared among different orders of potential with a fixed Pade approximant. Note that in Fig[2] (b), in the whole 
range of the figure, the Lo curve almost identically coincides with the Nnlo curve. 

Globally, the improvement with chiral orders is obvious: Nnlo prediction (solid line) is better than Nlo prediction 
(dashed line), and Nlo prediction is better than Lo (dotted line). There are also some interesting details: From these 
figures, we could see that in the higher energy region, all the Nlo curves have larger deviation from PWA data than 
the Nnlo curves. Some times, they are even worse than the Lo curves (C.f. Fig[2] (b) and (c)). 
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Here, we note that the solid curve in FigQ](a) seems puzzling. With only leading order potential (Vi ff ~ % Iq^+mF' 
plus a contact term V c = Co), one obtains pretty good predictions for the phase shifts, especially at the higher 
energies. The reason lies in the Pade approximant of Q, Re(£?) ~ (1 + 7^P 2 )/(f^ + ^-P 2 )); which is fixed 'physically' 
and effectively 'induces' higher order interactions upon iteration. This result agrees with the findings in Ref.[14j, 
where a nonperturbative of T-matrix is obtained using a simple potential (V = Vi^(Lo) + V(2)(in Eq.([5]). Nlo), and 
the prediction of the phase shifts is surprisingly good in a wider range of energy after the nonperturbative divergences 
is removed through fitting. Our analysis above provides a simple explanation of this surprise: The nonperturbative 
renormalization is properly treated! A closer look at the lower energy regions reveals that the higher order predictions 
still dominate the lower order ones (Cf. Fig|5](c)). Generally, the lower order predictions could 'win' at higher energies 
only by chance. 



UTILITIES OF PADE APPROXIMANT AND DISCUSSIONS 



Now we have seen that after the nonperturbative prescription dependence is properly resolved (here realized through 
low energy region fitting), the EFT approach facilitates 'physical' prediction for low energy NN scattering, at least in 
Sq channel. Since no specification of regularization and renormalization is needed, the Pade approximant of Q defined 
in Eq. ([1]) in fact provides a universal parametrization of the renormalization prescription dependence of the T-matrix 
in nonperturbative regimes. Note that both the potential and the Pade approximant of Q could be systematically 
extended to higher orders in EFT. 

Comparing with previous results, we find that the Lo prediction with ReOQ) ~ (1 + —p 2 )/{ — + —p 2 ) (C.f, Fig[3] 

p~l ^0 ^0 f 

(c)) differs significantly from that given in Ref. 4] and looks better. This nontrivial difference in predictions at leading 
order reflects the importance of nonperturbative renormalization. However, at higher chiral orders, especially at Nnlo, 
our results show no obvious differences in comparison with Ref.[3]. That means, including higher order interactions 
would lessen or tend to remove the nontrivial nonperturbative renormalization prescription dependence. This is a 
marvellous fact, since the fundamental requisite in EFT application is that renormalization prescription dependence 
should decrease as higher order interactions are included. Therefore, in spite of being an approximation approach, the 
procedures described above substantially proved or ascertained the rationality and applicability of EFT method in 
nuclear forces in a very general context. This is in sheer contrast to most known approaches where a renormalization 
prescription must be specified and hence the exploration of the prescription dependence is apparently limited. 

So far we determined the Pade parameters through fitting with the potential defined by EGMQ where the contact 
couplings were determined within a special cutoff scheme. In principle, we should determine the couplings in a way 
that is more prescription- independent: fitting through the combined space [C...] 1J [?..., /i], which should lead to a 
better way for determining the EFT couplings. Now the Pade approximant provides us a convenient approach to do 
so without really carrying out the formidable task of loop integrations and renormalization to all orders. We will 
perform the investigations along this line in the future. We believe other utilities could be derived from this analysis 
and the parametrization Eq.|T|). 
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From the above results, it is also obvious that it is fairly sufficient to employ Pade up to Re((?) ~ (1 + T^P 2 )/(r^ + 
l^p 2 )- For some channels, say P-wave, it is often sufficient to use Ke(Q) w + g(i)p 2 , which will be demonstrated 
in a separate report. 



SUMMARY 



In summary, we performed a Pade analysis on a compact factor of the T-matrix for NN scattering so that the 
nonperturbative prescription dependence and related effects could be conveniently explored. Such analysis suggests 
a useful theoretical instrument as well as a general and prescription-independent approach to test the efficiency and 
rationality of the application of EFT methods in nonperturbative regimes. Some related literature were also explained 
and discussed in favor of our analysis. 
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(a) (b) (c) 

FIG. 1: Predictions for 1 So phase shifts with three different Pade approximants: (1) Re(tJ) ~ gr \ (dotted line); (2) Ke(Q) « 
5(0) + 5(i)P 2 (dashed line); (3) Re(£) « (1 + ^p 2 )/(^ + ^P 2 ) (solid line). The circles denote the PWA data[23]. (a) for Lo; 
(b) for Nlo; (c) for Nnlo . 
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(a) (b) (c) 

FIG. 2: Predictions for 1 5o phase shifts at different orders of potential: Dotted line for Lo, dashed line for Nlo and solid 
line for Nnlo. The circles denote the PWA data[23]. (a) for Re(£) w g m ; (b) for Re(£?) ~ ff(o) + ff(i)P 2 ; (c) for Re(G) ^ 

(l + ^ 2 )/(^ + V)- 



